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The /3-decay of Li into Li and d is studied using simple halo wave func- 



tions of ^^Li. The sensitivity of the transition probability is elucidated on a 
^^ \ description of the halo part and on a choice of the potential between ^Li and 

d. 

^ ' Experiment with radioactive nuclear beams led to the discovery |I| of light nuclei with 

unusually large matter radii near the neutron drip-line. The nucleus ^^Li has received much 

•i-H . 

r\ \ attention as one of those typical systems which have neutron halo-like structure. The direct 

information on the extended structure was derived from the enhancement of the interaction 
cross sections. Momentum distributions of fragments have offered another useful way of 
probing the correlated structure of the halo part and testing model descriptions. Beta- 
delayed particle emission from the nuclei near the drip-line also provides us with a unique 
probe for the halo structure. The ^He j3 decay into a and d with the branching ratio, 
Pd=(7.6±0.6)-10^^, for Ed > 350keV is a good example and has recently been studied by 
several authors ||^-§[ . It was stressed in ref . , among others, that this decay is quenched by 
a cancellation between the internal and halo contributions to the Gamow-Teller(GT) matrix 
element and shows an extreme sensitivity to the halo description up to large distances. The 



purpose of this letter is to do an explanatory analysis on the ^^Li /3 decay into ^Li and d. 
This decay attracts much attention and is expected to give important information on the 
halo structure of ^^Li. 

Most of the (3 decay(~97%) of ^^Li takes place to low-lying states in ^^Be. A small branch 
of the P decay involves delayed n, t and a emission. Because of the experimental difficulty 
there is no direct evidence for the /^-delayed d emission from ^^Li . There is, however, 
only one experiment that has measured the t+{d) spectrum down to 350 keV although 
no separation was made between t and d. The experiment indicates that two separate 
components appear in the spectrum. The intensity at lower energy may be attributable to 
the branch of the deuterons because the Q values for ^Li+t and ^Li+d decay channels are 
4.96 MeV and 2.76 MeV, respectively. The lower limit on the branching ratio determined 
from the experiment is Pi+(rf) = (1.8± 0.3)-10^^ and a tentative value for the deuteron branch 
is obtained as Pd=(1.0± 0.3)-10-^. 

Despite numerous investigations the understanding of the structure of ^^Li is still imper- 
fect. This is in sharp contrast to the case of ^He in which a ^He+n+n three-body model 
works nicely. The three-body configuration of ^Li+ n+n is certainly an important ingredient 
for the ^^Li structure but there are at least two basic points to be settled before a more 
satisfactory description of ^^Li is made, i.e., the role of the polarizability of ^Li and the 
feature of the ^Li-n interaction. At present we have to use those very simple wave functions 
for ^^Li which assume a ^Li core, and our study is therefore of qualitative nature. 

We assume that the /3-delayed d emission proceeds through the so-called direct decay, 
i.e., a transition from ^^Li to ^Li+d continuum. In this (3 decay the isospin changes and the 
total angular momentum, Jj, in the final channel may be J/=l/2, 3/2, or 5/2. Since the 
decay Q-value is small, we assume the relative motion between ^Li and d is an s-wave and 
independent of Jj. Under these assumptions the decay transition probability per time and 
energy units, dW/dE, is calculated by g,^ 



where m is the electron mass, v the relative velocity between ^Li and d, (7/3=2.996x10"^^ 
the dimensionless /5-decay constant, and / is the Fermi integral. The GT reduced matrix 
element is expressed as 

BgAE) = 6A2 < SEltPeff >\ (2) 

where A = —1.25 is the ratio of the axial- vector to vector coupling constants. The radial 
part of the ^Li-c? relative motion function qe is normalized as 

gE{R) -^ R-\Fo{kR)cos6 + Go{kR)sm6), (3) 

where k is the wave number of the relative motion with energy E and 6 is the s-wave phase 
shift at energy E. The effective wave function ipeff is defined as 

/oo 
F,ir)xoir,Rydr, (4) 

u 

where Fd{r) is the radial part of the deuteron wave function. The function xo(^5-R) is the 
s-wave part in both of the coordinates, r and R, of a two-neutron reduced amplitude x(r, R) 
defined by 

x{r,R) = (55)1/2 ^ ^j^,^{%i),SMs = 00, TM^ = ll\^i^C'U) >, (5) 

where r is the relative coordinate between the two neutrons and R is the relative coordinate 
between ^Li and the center-of-mass of the two neutrons. 

We used the ^Li-d relative motion function qe generated from the folding potential of 
Watanabe type 



ViR) = 1 1 F,\r) 



U,{\R+'^\) + Un{\R-'^\ 



dr, (6) 



where Up{Un) is a p{n)-^Li optical potential of Woods-Saxon form. We assumed Up=Un and 
took into account only the real part of the central potential(Vo=53 MeV). Dependence of the 
transition probability on the optical potential parameters, the diffuseness a and the radius 
parameter tq, will be discussed later. The Coulomb potential of uniform charge distribution 



was used. As stated in the beginning our purpose is not to make a quantitative prediction 
but to learn the sensitivity of the spectrum and the branching ratio on the halo structure and 
the continuum wave function. We consider following simple two-neutron reduced amplitudes 
which were employed in the analysis of the momentum distribution of a ^Li fragment 
arising from the ^^Li reaction: (i) a (pi/2)j=o harmonic-oscillator shell model(SM) wave 
function, (ii) a di-neutron cluster model(CM) wave function, and (iii) a hybrid model(HM) 
wave function. For each case the function xo is expressed in terms of a nodeless harmonic 
oscillator function, F{b,r) = 2^-^^%-^/^ exp{-r^ /2b^) g: 

X^ir, R) = Ib^'ir' - 4R')F{2'/X, r)F{2-'/\, R), (7) 

X^^{r,R) = F{h,r)F{h,R). (8) 

The parameters, bo and (61,^2), are chosen to reproduce the measured interaction cross 
section of ^^Li on ^^C at an incident energy of 800 MeV/nucleon. See Table 1 of ref. 
1^ for the values. These two-neutron halo wave functions, however, cannot explain the 
observed momentum distribution of ^Li in the (^^Li, ^Li) reaction. In case (iii) a simple 
wave function that reproduces the momentum distribution was constructed [^] by taking 
a linear combination of the shell model and cluster model wave functions: x^*^(r, i?) = 
eiXo^'^ {r, R) + ^2X0^^ {^^ R) with ei = 0.79 and 62 = —0.41. The hybrid model is considered 
"best" among three cases. 

Fig.l illustrates how the functions qe and ipeff contribute to the GT matrix element for 
E=0.5 MeV ( a=0.6 fm, ro=1.2 fm). The integral, I{R) = Jq gE{R') i^effiR') R'^dR', for 
three effective wave functions is displayed in the lower part of the figure. The limit of I(R) 
for R ^ 00 appears in eq.(2). The GT value gains most of the contribution in the region of 
i?=4~10 fm. It has negligible contribution from the interior part of the wave function, which 
is understood by the fact that the ^Li-d relative motion function exhibits more oscillatory 
behavior in the interaction region than the a-d case. We can conclude that the (3 decay 
of ^^Li into ^Li and d probes mainly the halo part of the ^^Li wave function. This is in 



sharp contrast to the case of ^He, where the cancellation in the GT matrix element occurred 
between the internal and external parts of the ^He wave function. 

We have compared the energy dependence of the transition probability for three model 
wave functions. All three cases yielded similar spectrum although the magnitudes of the 
transition probability differed within a factor of ten at the peak values. 

Table | lists the branching ratio calculated with the HM wave function for various values 
of the potential parameters, a and tq. Since no ^Li+d scattering data are available, we cannot 
select a good set of the parameters. The potential with ro=1.3 fm is ruled out, however, 
because it leads to a smaller branching ratio than expected from experiment. Increasing 
To or a leads to a stronger potential. The branching ratio of order 10~^ is obtained when 
the third node of IqeI is pulled inside at shorter distance and the overlap integral with ip^ff 
becomes smaller as the potential becomes more attractive. Fig. 2 shows the dependence of 
the transition probabilities on the potential parameters. The energy at which the peak of 
the spectrum occurs is dependent on the choice of the potential. Measurements with good 
statistics down to the lower energy will be very informative in this respect. 

In conclusion we have examined the /3 decay of ^^Li into ^Li and d in order to shed light 
on the role of the halo structure of ^^Li. We find that the /9-decay matrix elements are 
determined practically by the halo part of distances, R > 4 fm. The branching ratio and 
the ^Li-d relative motion energy at which the transition probability reaches a maximum are 
sensitive to the choice of the optical potential. A standard set of the potential yields the 
branching ratio of order 10~^. 
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TABLES 
TABLE L Dependence of the branching ratio of the ^^Li (3 decay into ^Li and d on a set of the 
parameters, a and ro (in fm ), of the optical potential. The hybrid model wave function is used. 
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FIGURES 
FIG. 1. The contribution of the initial and final wave functions to the GT matrix element. 

Absolute values of RgE^R) for E = 0.5 MeV and of Ripeff{R) and the integral I{R) are drawn. 

See text for the definition of I{R)- 

FIG. 2. The transition probability per time and energy units as a function of the relative 
motion energy E for different sets of the optical potential parameters. The values of a and vq are 
given in units of fm. The hybrid model wave function is used. 



